Independent particle descriptions of tunneling from a many-body perspective 
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Currents across thin insulators are commonly taken as single electrons moving across classically 
forbidden regions; this independent particle picture is well-known to describe most tunneling phe- 
nomena. Examining quantum transport from a different perspective, i.e., by explicit treatment of 
electron-electron interactions, we evaluate different single particle approximations with specific ap- 
plication to tunneling in metal-molecule-metal junctions. We find maximizing the overlap of a Slater 
determinant composed of single particle states to the many-body current-carrying state is more im- 
portant than energy minimization for defining single particle approximations in a system with open 
boundary conditions. Thus the most suitable single particle effective potential is not one commonly 
in use by electronic structure methods, such as the Hartree-Fock or Kohn-Sham approximations. 

PACS numbers: 05.60.Gg 73.40.Rw 73.63.-b 31.25.-v 



Describing quantum transport from first-principles has 
proven to be a challenging task; debate continues as to 
the proper theoretical approach for treating electron cur- 
rents across metal-molecule-metal junctions (MMJs) 0, 
Hi IE Q H) IE Q' Tunneling is an archetype of quan- 
tum behavior with implications for all branches of mod- 
ern physics, as well as chemical and biological processes 
dominated by electron transfer. The ability to accu- 
rately describe tunneling currents on the nanoscale is 
important for both fundamental science and technology 
perspectives. The discussion focusses on the suitability 
of various Kohn-Sham density functional theory imple- 
mentations (KS-DFT) US 110, or possibly Hartree- 
Fock method 0, @i ^o apply in combination with the 
one-body non-equilibrium Green's function (NEGF) ap- 
proach 0, |E|. In fact though, there is no estabhshed 
criterion for selecting a single particle Hamiltonian to be 
used in transport calculations. 

On the other hand, evidence has been mounting that 
electron transport can be sensitive to many-body ex- 
change and correlation 0,0,0- For example, applying 
NEGF together with static KS-DFT leads to overestima- 
tion of the current in MMJs, in some cases, up to orders 
of magnitude when compared to experiments [j| . Several 
reasons for discrepancies have been suggested: DFT's un- 
derestimation of occupied/unoccupied state separation 
results in too high currents 0, , conversely the over- 
estimation of the gap in the HF approximation yields 
currents too low !3| ; DFT exchange-correlation function- 
als do not accurately reflect potential profiles [lol : ot her 
methods such as time-dependent current DFT 'll' or 
density matrix functional theory are simply better suited 
for the treatment of electronic currents. Irrespective of 
the proposals to explain/remedy single particle tunneling 
descriptions, it is necessary to understand the physics de- 
riving from a genuine many-body formalism to identify 
the source of the discrepancies. 

In this Rapid Communication we investigate the ex- 



tent of correlations beyond the single particle picture 
and identify conditions for defining a "best" independent 
particle model, using a recently formulated many-body 
quantum transport approach [5| . We demonstrate a sin- 
gle electron description of tunneling deriving from maxi- 
mizing the overlap of a single Slater determinant with the 
true many-body current carrying state, and show that it 
remains a good approximation even after conventional 
independent particle models fail. 

Most theoretical studies of quantum transport begin 
with the use of a Slater determinant of single particle 
states to model a many-body current-carrying wavefunc- 
tion Ti'l. Our recent work |3| finds this picture becomes 
markedly less valid near or above a resonance. Here, we 
concentrate on non-resonant tunnel models. Independent 
particle models are appealing because as the wavefunc- 
tion can be represented as a single Slater determinant, 
the resulting physical model can readily be pictured and 
computed '4|ll3|. Their limitation is neglect of electron 
correlations; below we introduce a correlation measure to 
quantify when an independent particle model holds. 

For our study, MMJs are an ideal test case as the 
molecule acts like an insulator when there is little charge 
transfer or hybridization with the electrodes, there is a 
large set of experimental observations to validate the cal- 
culations, and MMJs may be modelled with a relatively 
small number of atoms allowing their electronic struc- 
ture to be accurately calculated with many-body meth- 
ods. We choose alkane chains (C2H4)„ (inset of Fig. 
and for comparative purposes also examine silicon hy- 
dride molecules (Si2H4)„ (silanes), as replacing the car- 
bon atoms with silicon reduces the occupied/unoccupied 
state separation and increases the degree of correlation. 

Details of the transport method are given in ref. |E A 
constraint Ansatz using the Wigner function is made to 
incorporate open-system boundary conditions for calcu- 
lation of the reduced density matrix on a device region 
(typically the molecule plus part of the electrodes) at sev- 
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number of carbons (2n) 



FIG. 1: Tunnel resistance-increase exponential law. Expla- 
nation of the standard basis set notations is given in the text. 
Inset: Typical geometry of the studied Autip-S-(C2ll4)„-S- 
Autip molecular junctions (n = 3). 

eral values of applied voltage. The procedure results in 
the best approximation to the density matrix on a region 
subject to reproducing known system observables in ac- 
cord with the principle of maximum entropy at the zero 
temperature limit. The density matrix is calculated from 
eigenfunctions of a constrained many-electron Hamilto- 
nian. A strength of our approach is that it allows for the 
expansion of many-body states in terms of a complete set 
of configurations 

\^) = co|vI/o)+^c^|^n+EE4'l^S'') + ---(l) 

i,a i<j a<b 

|4'o) refers to a reference state composed of TV low- 
est single particle states, \'^f) d^'fj') . . .) denotes singly- 
(doubly- . . .) excited configurations generated by substi- 
tuting the i-th (j-th, . . .) occupied single particle state 
with the a-th (6-th, . . .) single particle excitation. In- 
dexed l^*) are spin-projected Slater determinants or con- 
figuration state functions (CSFs), and |^'o) is the ground- 
state HF determinant in the absence of an applied elec- 
tric field. The symmetry group of our contact-molecule- 
contact subsystem is C2h (inset of Fig. P) whose ground 
state eigenfunction at y = is a singlet with Ag sym- 
metry. The bias field polarizes the molecular subsys- 
tem and mixes states of singlet Bu symmetry via dipole 
terms. Due to the combinatorial nature of generating 
possible configurations, there is a computationally pro- 
hibitive number of excited state determinants. To cal- 
culate the many-body tunneling wavefunction |v1/mb) we 
approximate the full expansion by taking the most 
significant CSFs, as computed from a Monte Carlo sam- 
pling of the many-electron expansion space |l4l | with the 
addition of all singly excited CSFs. Notably, this ap- 
proximation is not excitation limited and allows for a 
high degree of electron correlation . 

In contrast to the many-body wavefunction, a single 
Slater determinant I'I'i-dct) is sought in independent par- 
ticle models. If this is initially built from a set of unper- 
turbed single particle states, e.g., the Hartree-Fock \'i'o) 



as above, and a perturbation such as an electric field 
is introduced, the resulting many-body expansion of the 
perturbed single determinant state has the special fea- 
ture that the coefficients of the doubly- (triply-,...) ex- 
cited configurations will be appropriate products of the 
coefficients of the singly-excited configurations 

That this is the most general determinant not orthogo- 
nal to l^'o) is seen by expanding the Thouless expression 

|*i-dct) = <^^P{J2i,a'^i^l^t)\'^o), where ajia^) creates 
(annihilates) the i-th orbital. From the many-body cal- 
culation we find that there is a broad voltage range where 
linear response is valid; then by taking terms up to second 
order in the cf in ((21 and applying the same Wigner func- 
tion transport scheme to this restricted form of the wave- 
function, we deduce the HF results to linear response in 
the bias. This permits comparison of the exact-exchange 
(but uncorrelated) and correlated descriptions of tunnel- 
ing within the same formalism. 

In Fig.^both correlated and uncorrelated tunnel resis- 
tance of the alkane-based molecular junctions are shown. 
The qualitative features are the exponential suppression 
R = Rqc^^" of the current with respect to the number 
of carbon atoms iVc = 2n (Fig. ^ and linear scaling with 
applied voltage. As observed in numerous experiments, 
both results are typical of a metal-insulator-metal sys- 
tem. Contact resistances i?o smd tunneling parameters (3 
calculated by least-square linear fits are relatively close 
in value with HF i?o and (3 deviating from the many- 
body results by 32% and 14%, respectively. The over- 
estimation of (3 by the HF approximation can be inter- 
preted in terms of the incorrect alignment of the virtual 
orbitals which leads to a larger highest occupied-lowest 
unoccupied molecular orbital gap. In contrast, the con- 
tact resistance i?o is underestimated as indicated by the 
zero crossing in Fig. ^ suggesting a stronger molecule- 
electrode coupling within HF. Similar behavior is found 
for the silane chains, with j3 ~ 0.18 per silicon and a 
contact resistance of i?o — OOOKfi from the many-body 
calculation. As expected, these molecular junctions are 
much more conducting yielding a value of /3 « 0.09A~^ 
compared to (3 k, 0.39A^^ for alkanes. But, as will be 
discussed in detail, the silanes display a higher degree of 
electron correlation and the HF deviations with respect 
to the many-body results are substantially larger than 
for the alkanes, increasing to 75% and 117% for Rq and 
(3, respectively. 

Quantitatively, the magnitudes of the low-bias resis- 
tance for each alkane molecule compares well with those 
reported in ref. apart from decane which appears to 
be experimentally poorly resolved, as the measured re- 
sult for n = 5 does not fit the extrapolated tunneling 
behaviour. Our estimated inverse decay length [3 falls in 
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Diamonds: Au -S-(C,H J -S-Au 

Hp ^ 2 4'ii tip 
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molecular length (2n) 

FIG. 2: Many-body correlations expressed via the one-body 
density matrix. These are vohage-independent with a stan- 
dard deviation much smaller than the symbol size. 



the lower part of the experimentally observed range of 
values scattered between around 0.5 and 1.0 per carbon 
atom III particular, Haiss et al. have reported 
P = 0.52 ± 0.05 per carbon but with higher resistance 
values (compare to values indicated in Fig.QJ. 

We can state that HF performs relatively well com- 
pared to a full many-body calculation for non-resonant 
transport in the alkanes, but eventually diverges from the 
many-body result as the molecular length increases. Due 
to the higher degree of correlation, HF is a poorer ap- 
proximation for the silanes, and also diverges with chain 
length as compared to the many-body calculation. In 
fig-E the HF results for the smaller but commonly used 
atomic orbital split valence with polarization (SV(P)) 
basis set for the carbon atoms is included (the larger 
aug-cc-pvDZ denotes a valence double zeta correlation- 
consistent basis set with polarization, augmented with 
diffuse functions). Within this approximation, the molec- 
ular gap becomes larger yielding a small increase in P and 
higher resistance. This result does not compare as well 
to the many-body calculation and casts doubts on results 
whose convergence has not been tested with respect to 
the completeness of the single particle basis. 

We have examined the importance of electron corre- 
lation by comparing separate calculations at the many- 
body and HF levels, however, the level of correlation may 
be directly quantified from the many-body calculation 
independently of HF results by examining the one-body 
density matrix 71. Lack of correlations is expressed by 
the condition that 71 may be derived from a single deter- 
minant if and only if it is idempotent, Tr[7^ — 71] = 0. 
For closed shell systems such as studied here, the degree 
of correlation C > or equivalently the deviation from 
idempotency - for the spin-traced density matrix pi may 
be defined in terms of its eigenvalues or natural occupa- 
tion numbers n,- as 
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molecular length (2n) 

FIG. 3: (Left axis) Contribution of correlations to the current 
measured via Eq. |3 (Right axis) Current percentage not car- 
ried by the leading determinant defined as the A'^-particle con- 
figuration with highest natural occupancies. Residual noise 
voltage-dependence is indicated by error bars. 



As < rii < 2, an upper bound can be deduced: 
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Here N is the total number of electrons in the device 
region and Nt is the size of the single particle basis. The 
magnitude of C is a direct measure of the deviation of the 
actual one-body density matrix from an uncorrelated one 
deriving from a determinantal wavefunction; C = for a 
single determinant and C = Cmax when nt = N /Nt- 

In fig. 12 we plot the correlation measure C scaled to 
its upper limit as a function of chain length for I^'mb) 
(note that it vanishes for |4'i_dot))- Remarkably, as a 
function of voltage (not shown) the correlation measure 
is constant and close to the uncorrelated limit for any 
given oligomer. In fact, C monotonically reduces with in- 
creasing length possibly approaching a lower limit set by 
an infinite polymer for both the alkanes and the silanes; 
the latter system displaying more electronic correlation. 
Hence, in contrast to the HF predictions, the density 
matrix reveals there exists a single Slater determinant 
description that is equally valid for all chain lengths. 

A formal measure of correlation for any expecta- 
tion value (v&MBl-f I^mb) can be devised by using the 
fact that the spin-traced two-body density matrix p2 
factorizes in terms of pi only in the one-determinant 
(uncorrelated) approximation. We write p2{ij;kl) = 
(1/2) (pi {ik)pi (jl) - pi {il)pi ijk)/2) + X{ij; kl), where the 
first term is familiar from HF theory and A denotes any 
deviations when pi is not idempotent. Here i and k are 
indices of the first electron, j and I that of the second 
and ^ y02(u; u) = N{N — l)/2. Tracing over the coor- 
dinates of the second electron, we deduce a partition of 
pi — 1/(A^— l)[7V/9i— pi xpi/2]-|-A into uncorrelated and 
correlated pieces. This gives us a decomposition of the 
current (/) = Tr[yOi/] flowing in the many-body wave- 
function as ImB = -^l-dct + Icorr, whcrC Icorr iucludcS 

the contributions from A and does not vanish when pi 
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voltage (V) 

FIG. 4: Comparison of the (a) Many-Body, (b) Hartree-Fock, 
and (c) leading natural orbital determinant current. 

derives from a correlated wavefunction. It follows 

^cor. = ^^Tr[(p2/2-pi)/]. (5) 

The percentage of current carried by the correlation 
terms is indicated in fig.|31by the filled symbols. In line 
with our previous observations, it is characterised by a 
single value for each oligomer. The correlation contribu- 
tions remains a small fraction of the total current as a 
function of applied voltage. These findings imply the ex- 
istence of an approximate single-determinantal current- 
carrying wavefunction that captures most of the contri- 
butions to electron flow even when the HF approximation 
begins to fail. 

The HF approximation determines the best orbitals 
for a l^'i-dct) that minimizes the energy. Since the cur- 
rent largely reflects the form of the wavefunction, it is 
plausible that an approximate state that maximizes over- 
lap with the current carrying wavefunction may improve 
the orbital description of an independent particle model. 
Truncated configuration expansions built from the eigen- 
functions of the one-body density matrix, the natural or- 
bitals (NOs), are known to maximize overlap to 
Since there are only 7V/2 NOs with large occupation num- 
bers rii — in agreement with the calculated small corre- 
lation measure C in our MMJs — we use the single deter- 
minant made by doubly occupying the orbitals with occu- 
pation numbers close to 2; note that an exact treatment 
of maximizing the overlap for |vE'i_dct) yields the Brueck- 
ner determinant. We find that the current carried /no 
is close to /mb as shown in fig. |31 where the percentage 
of the remaining differences is plotted versus molecular 
length. Remarkably, these deviations are of the same or- 
der as the correlation contributions determined by eq. |S1 
Indeed, in fig.^jthe current- voltage characteristics clearly 
demonstrate that the "maximum overlap model" yields 
almost indistinguishable results from the many-body cal- 
culations independent of residual electron correlations; 
the HF single-electron picture already shows significant 
deviation for pentane (fig.EIa)) which grows much larger 



for the silanes (fig.Efb)). The use of the natural orbitals 
in defining |4'i_dGt) acquires a practical significance in 
view of recent advances allowing their construction from 
one-electron equations with an effective potential [iflj l . 

In summary, we have shown within a many-body 
scheme how an independent particle description emerges 
for electrons tunneling across the barrier of insulating 
materials, due to zero-bias correlations being such that 
a single determinant wavefunction is appropriate. Our 
main conclusion is that a strong selection criterion for a 
single particle transport model is to maximize the over- 
lap with the many-body state; this does not yield the 
Hartree-Fock determinant. Indeed, this Slater determi- 
nant outperforms HF substantially. As most KS-DFT 
implementations to date have recently been shown to be 
performed at essentially the Hartree level the indepen- 
dent particle picture we have established here is better 
than both present models and, hence, valid beyond the 
limits of conventional electronic structure methods for 
quantum transport. 
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